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We wish to consider in this report the large N limit of a particular matrix model introduced by 
Myers describing D-brane physics in the presence of an RR flux background. At finite N, fuzzy 
spheres appear naturally as non-trivial solutions to this matrix model and have been extensively 
\l ' studied. In this report, we wish to demonstrate several new classes of solutions which appear in 

^^ ' the large N limit, corresponding to the fuzzy cylinder,the fuzzy plane and a warped fuzzy plane. 

The latter two solutions arise from a possible "central extension" to our model that arises after we 
account for non-trivial issues involved in the large N limit. As is the case for finite N, these new 
O ' solutions are to be interpreted as constituent DO-branes forming D2 bound states describing new 

O , fuzzy geometries. 

! I. INTRODUCTION 

(N 

^\ • Non commutative geometry, not withstanding its intrinsic interest T] ^ I3J has attracted much attention recently 
T— I ' through its apearance in string theory and M-thcory (sec 4J 5] 6] 7] and references therein). One of the more 
04 familiar manifestations of non commutative geometry in this context arises when one considers the dynamics of 

^D ' D-branes [g in various backgrounds. There is a natural matrix model description of this dynamics, and in general 
one expects the solutions of this matrix model to correspond to various fuzzy geometries. In our report, we wish to 
study a particular matrix model that arises when we consider the dynamics of DO branes in the presence of an RR 
flux background 0E3- Although this model arises in various other contexts (see 01)j we will only be interested in 
this model in so far as it provides for a framework in which to study the dynamics of noncommutative spaces. The 
I \ interpretation we prefer is that which arises from Myers' 5] derivation of this matrix model: that the dynamics and 
jr~^' energetics of this model correspond to various bound DO brane configurations. The dimension of the matrices we 
(-H I consider corresponds to the number of DO branes we have in our system, and so the large N limit corresponds to 
• • . taking the limit of a large number of DO branes. This is certainly a reasonable scenario to consider from the 

perspective of early universe physics. Indeed it is with with one eye on spacetime physics that we undertake this 
K^ , study, as we wish to use this report as a point of departure for future work Q , where there are preliminary 
\^ • indications of dynamical D-brane topology change in this large N limit. 

a 

Although it might seem at first sight that the large N limit involves a loss of calculational ease when considering the 
dynamics of the system, we find that after careful consideration the situation is not that much different from the 
case at finite N. The few important diffrences we will uncover only serve to enrich the problem, and allow us 
solutions that we would not have had at finite N. In particular, we find that such considerations yield solutions that 
correspond to the DO-branes configured as non-commuting planes, and non-commuting cylinders. The possibility of 
transitions between these new geometries will be taken up in a follow up report 13 • 
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II. THE MODEL: A FIRST PASS 

The matrix model we will be studying is described by the following action (QIhIQIH): 

5 = To fdtTr(^^X,X, + ^[X,, Xj][X,,X,] - '-Ke,jkX,[X,,Xk]) (1) 

Following the conventions of j4|, we work in units where 27rQ;' = 1. The Xi are taken to be TV x A^ traceless, 
hermitian matrices and Tq = \f^ j gs is the DO brane tension. The index % runs from 1 to 3. In addition to this 
action, we have the Gauss law constraint: 

\X,,X,\=Q (2) 

Which arises from the Aq equation of motion for the RR gauge field. Note that the last term in the action is a 
Chern-Simons term which induces interactions between the DO-branes through the 4-form flux, which assumes the 
vacuum expectation value ^0123 = — SKe^fc. This term was deduced by Myers by demanding consistency of the 
D-brane action with T-duality |5(. 

Let us for the moment forget dynamics, and concentrate on the static solutions of this model. Varying the action 
gives us the static equation of motion: 

[X,,([X„X,]-^«ey■feXfc)]=0 (3) 

From which we can immediately deduce two classes of solutions: commuting matrices ([X^, X.j\ ~ 0) and fuzzy two 
spheres- 

\X,,X.^=iKe^kXk (4) 

Any representation of the lie algebra of SU(2), irreducible or otherwise, given by the generators { Ji} furnishes the 
solution (see [l3| for a review of the non commutative 2-sphere) : 

Xi = nJi (5) 

Where for an irrep of su(2), the casimir operator XiXi = K^j{j + 1) gives us the radius squared of the 2-sphere and 
the dimension of the irrep is given by A^ = 2j + 1. Clearly, a reducible representation can be expressed as a direct 
sum of irrcps: 

X, = ®,K jf ) (6) 

The summation over r runs over the irreps included in this representation and describes a solution of multiply 
superimposed fuzzy 2-spheres. As before, the radius of each constituent fuzzy 2-sphere is given by R^ = K^jrijr + 1) 
and the dimension of the representation is now given N = '^j,{'2jr + 1). 
The energy of these solutions is given by their potential energy: 

V = Tr(^-^[X,,Xj][X,,X,] + ^Ke,,kX,[X,,Xk]) (7) 

From which we see that commuting matrices correspond to zero energy solutions. Without the presence of the 
Chern-Simons term in our model, these solutions would span the lowest energy configurations. However the addition 
of the Chern-Simons term modifies things drastically by permitting fuzzy sphere solutions, which have a negative 
energy: 

E = -Tok^IY.-^-^^^ + 1)(2> + 1) (8) 



So we see that the fuzzy two spheres describe a bound state of DO branes, with the irreducible representations 
corresponding to the lowest energy configuration for a given N. The observation that reducible representations have 
greater energy than irreducible ones has motivated the belief that this model could describe topology changing 
physics through transitions between these states. The simplest example of a transition that this model might be 
capable of would be when two separate spheres meld into one. Observations based on studying the energetics of this 
model iJllalH have yielded promising preliminary results, but an explicit solution interpolating between two 
topologies is still lacking. We wish to use the results collected in this paper to further investigate the existence of 
topological transtions in a future report 9], and in the process hopefuly highlight why these solutions cannot be 
taken for granted even if the energetics of the model strongly imply the possibility of topology change. 

Returning to the static equations of motion (jS)) , it might appear that we have exhausted all the possible static 
solutions, however a slight recasting of the problem will show that this is not the case. Performing the change of 
variable: 

X+ = Xi + 1X2 , X_^Xi~ 1X2 ; xl^X. (9) 

The equations of motion become: 

X+ = ^[X+,[X+,X_]] + [X,,[X+,Xs] + 2kX+] (10) 

X- = i[X_,[X_,X+]] + [X3,[X^,X3]-2kX^] (11) 

^3 = i[X+,[X3,X_]] + ^[X_,[Xs,X+]] + k[X+,X_] (12) 

And the Gauss law condition becomes: 

[^,^3] + i[X+,X_] + i[X+,X_] = (13) 

Since pO|l and pi|l arc adjoints of each other, we are left with only two independent matrix equations to satisfy. 
The fact that one of the equations governs the motion of a non-hermitian matrix, which in general has a hermitian 
and an anti-hermitian part, accounts for the two hermitian equations we had previously. We now make the ansatz: 

[X+, X_] = AX3 , [X3, X±] = ±9X± (14) 

For which eqs. (^HJ-lEJ imply 

X+ = {2K-e -\/2)0X+ (15) 

X3 = (K-e)AX3 (16) 

Thus in the static case, we can immediately read off three distinct solutions. A = = is one solution which 
corresponds to commuting matrices. A = 2k, 9 = k is another, which corresponds to the already known fuzzy two 
sphere solutions. It should be clear that the change of variable we made is identical to constructing the usual raising 
and lowering operators out of Xi and X2 , and the algebra: 

[X+, X_] - 2^X3 , [X3, X±] = ±kX± (17) 

is none other than that of su(2). The third solution is given by A = 0, 6* = 2k: 

[X+,X4 = , [X3,X±] = ±2kX± (18) 



Which corresponds to the fuzzy cyhnder MTl provided we impose the additional constraint X+X- = p^I. Since the 
non-commutative cyhnder is not as famihar a fuzzy geometry as the sphere or the plane, we briefly outline how it is 
constructed, following the treatment of [Ij. We wish to interject at this point that alarm bells should already be 
ringing, as l|18|) describes the lie algebra of a non semi-simple group (re-expressing the above in terms of Xi, X2, X3 
we see that Xi and X2 form a proper ideal). Hence we will neccesarily be dealing with infinite dimensional 
representations if we want to preserve the hcrmitian nature of the generators (i.e. if we want to work in a unitary 
representation)- which is required by the very nature of our model. It is not immediately clear what equations of 
motion the action (^) will yield as N tends to infinity, as crucial in our derivation of the equations of motion was the 
cyclic property of the trace, which allowed us to bring all of the variations of the Xi to one side of the trace |2J|. The 
cyclic property of the trace cannot be taken forgranted in infinite dimensions, and for unbounded operators it 
generally fails. We shall see in the next section that rather remarkably, we can be forgiven for proceeding as we are 
at the minute. Indeed, there we will discover a possible fourth class of solutions after accounting for the subtleties of 
the infinite N limit. However for the present purposes we will take these issues forgranted and justify our treatment 
here later. 

Returning to the fuzzy cylinder, it helps if we begin with the commutative cylinder, described by a non compact 
coordinate t e R and an angular coordinate e [0, 27r]. We see that functions on the cylinder are spanned by the 
basis: 

{e^^nez ; /(r, 0) - ^ c„(r)e'"^ (19) 

n 

Defining x+ — pe**^, X- — pe"*"^, where p is the radius of the cylinder, we can generate all fmictions on the cylinder 
by a power series in these variables. Furthermore, we have the relation x+x^ — p^. The Poisson structure of the 
cylinder is defined by the brackets: 

dfdg dgdf 

where / and g are arbitrary functions on the cylinder. By the Leibniz rule, all Poisson brackets can be generated 
from the following elementary brackets: 

{t,x±\ = ±ix± , {x+.x^} = (21) 

The relationship x+X- — p^ is preserved by these brackets (i.e. it is central): 

{t,x+X-} — {x±,x+x^} = (22) 

Derivatives of functions on the cylinder can now be effected by the action of the brackets: 

dll = {r, {r, /}} ; d^f = -l{a;_, {x+, /}} (23) 

To obtain the non-commutative cylinder, one 'quantizes' the Poisson structure through the prescription {, } — > 4[, ] 
where A is the non-commutativity parameter. Although one feels that this procedure shouldn't raise any suspicion 
due to its similarity with canonical quantization, it is in fact a very well defined and well motivated prescription for 
constructing fuzzy geometries ^5]. In fact a quick and easy way of seeing how this is plausible, is to realise that 
because of the remarkable formula [A, BC] — [A. B]C + B[A, C], the above quantization prescription reproduces all 
Poisson brackets of functions, and hence their derivatives, in the form of commutator brackets acting on the 
corresponding operator functions, up to ordering. Returning to the problem at hand, we can thus write down the 
structure relations of the fuzzy cylinder as: 

[t, a;±] ~ ±Ax± , [a;+,x_] — , x+x- = p (24) 

Thus if we take A = 2k, we see that the solution to H18|l is indeed the fuzzy cylinder if we identify X3 with r and X± 
with x± . From the algebra H18|l one can immediately glean that the spectrum of X^ is integer multiples of 2k and 



the X± act as raising and lowering operators in the basis where X^ is diagonal. In fact because of the relation 
X^X_ — p^, we can conclude that the action of X^ is that of p times the elementary shift operator: 

Xsln) = 2Kn\n) -^ X+\n) = p\n + 1) (25) 

Where it should be clear that X_ performs the opposite shift. It is this qualitative similarity with the situation for 
the 2-sphere (where of course the raising operator has a more complicated dependence on the state it is acting on 
through the Clebsch-Gordon coefficients) that motivates us to investigate in Q, the possibility of a transition 
between these two fuzzy geometries. 

Thus in addition to forming bound states corresponding to fuzzy 2-spheres, a large number of DO branes in the 
presence of RR flux can also form a bound state corresponding to the fuzzy cylinder. However there is more to this 
model yet, in the next section we will confront the issues involved in dealing with infinite dimensional matrix 
actions, and find through our efforts yet another distinct fuzzy geometry admitted by this matrix model- the 
noncommutative plane. 

III. INFINITE DIMENSIONAL MATRIX DYNAMICS 

In deriving the equations of motion from the action (^ defining this matrix model, we had to make use of two 
crucial aspects of the trace. The first aspect is that it is cyclic- for an arbitrary variation SXi will not commute with 
the Xi themselves, hence we need this cyclicity to bring all variations to one side of the trace. The second aspect we 
need is the positive definiteness of the trace norm, which allows us to deduce the equations of motion from the 
statement that the action is extremized by all variations in the fields: 

Tr(^-^SX?j ^OySX, ^ ^ = (26) 

Now in the case where the matrices we are dealing with are infinite dimensional, and the trace is over a suitable 
Hilbert space, we cannot take the cyclic property of the trace forgranted. Certainly the trace remains cyclic for 
bounded operators, but when we are dealing with unbounded operators this is not the case (e.g. X, where the 
spectrum is the real line). The most famous example comes from the canonical pair X and P with [X, P] = i. It 
should be apparent that if were we to work in a basis where any one of the pair is diagonal, then clearly 
TrXP ^ TrPX . In fact far from being a mysterious mathematical oddity, the origin of the non-cyclicity of the 
trace lies in the mundane fact that one cannot interchange the order of integration over a domain where the function 
we are integrating is unbounded. 

This immediately begs the question of how we should proceed deriving the equations of motion if we want to drop 
the assumption that the matrices we are working with are finite dimensional. The answer is just as immediate and 
for the most part, the end result is not that different in many respects except for one, which will facilitate a "central 
extension" to our matrix model. Before we continue however, we wish to point out that the only regularization 
procedure we'd ever need to carry out in order to make the transition to infinite dimensional matrices, is to 
incorporate a normalization factor into the trace such that TrI — 1. This is only ncccesary when studying the 
energetics of the system, and has no physical significance in terms of a cut-off length scale. Rather, all this serves to 
do is to factor out the divergent behaviour of the energy such that we can meaningfully compare the energetics of 
different configurations through their energy densities, which always remains finite. 

Varying the action Q as it stands, we end up with the following: 

dS = ToJdtTr(^^6X,X, + ^X,SX, + ^[6X„Xj][X,,Xj] + ^[X,,Xj][SX„Xj] (27) 

— -iriijkoXi[Xj,Xk\ —eijkXi[5Xj,Xk\\ (28) 

Where for finite dimensional matrices, we compensated for the non-commutativity of the fields and their variations 
by using the cyclic property of the trace. In the present case, we proceed as follows. We assume that the matrices 



are now operators over some separable Hilbert space Ti. The basis set of this space \q),q e R satisfies the following 
properties: 

j dq\q){q\=I- {q\q') = 6{q - q') (29) 

From which we deduce that the matrix elements: 

0{u,v) -.^ {u\0\v) (30) 

define completely an operator over this space. The trace of an operator is given by: 



dq{q\0\q) = / dqO{q,q) (31) 

And the product of two operators is given by: 



(Oi02)(p,g) = {p\OP\q)^ J drOi{p,r)02{r,q) (32) 

Where hermitian operator matrix elements satisfy: 

0^p,q) = 0*{q,p) (33) 

One can see from the trace of the commutator of a pair of operators: 

Tr[0,,02]= J dq{q\[Oi,02]\q) ^ J dqdp(Oi{q,p)02{p,q)-02{q,p)Oi{p,q)) (34) 

that although the integration variables are dummy variables (and hence can be interchanged through relabeling), 
the order of inegration cannot be interchaned. If this were the case, the above expression would vanish identically. 
Hence in all that follows, we will meticulously preserve the order of integration when writing out operator products, 
with the order of integration is to be read from right to left. 

Before we can proceed further, we have to address two important issues that constrain this model. The first issue is 
that of the Gauss law constraint (01, and the second is that our matrices are constrained to be hermitian, and in the 
event that the translational symmetry of this model is broken, of constant trace (we demonstrate further on that 
this symmetry is spontaneously broken in the infinite dimensional limit). We will approach the first issue in the 
usual manner of introducing auxilliary fields. Although one might be tempted to tackle the second issue in a similar 
manner, thus permitting arbitrary field variations when applying the variational principle, we choose to account for 
it more directly by only allowing variations that are consistent with preserving the hermitian and/or constant trace 
nature of the matrices. This not only does away with what would be a proliferation of auxilliary fields, but turns out 
to be rather easy to implement. 

Reconsider the constraint: 

Because of the particular form of our action, we have the luxury of two different choices in implimenting this 
constraint. One could either introduce an auxilliary matrix, A such that we add the following to the action: 

dtTr{X[X,,X,]) (35) 



Where requiring the action to be stationary under arbitrary variations of the elements of A give us the desired 
constraint. We can alternatively introduce an auxilliary function of time, A and add the following term to the action: 

dt\Tr{[X„X,][X,,X,]) (36) 

Requiring the action to be stationary under variations of A(i) implies: 

-Tr{[X^,X^][Xj,Xj]) = (37) 

Since this equation is of the form TrO^O, by the positivity of the trace norm we can again conclude the desired 
constraint. It turns out that the latter approach is more convenient as it avoids having to solve for the constraint 
explicitly. 

Turning now to the second issue constraining our model, we see that if our variations are to preserve the hermeticity 
of the matrices, then the variations themselves have to be infinitessimal hermitian operators. In terms of matrix 
elements in our Hilbert space this implies: 

SX,{q,p)^SX*{p,q) (38) 

An arbitrary variation that satisfies this constraint is then given by linear combinations of the following: 

Si,e,u.vXi{q,p) = Sii{eS{q - u)S{p - v) + e*S{q - v)S{p - u)} (39) 

Where e is an arbitrary phase. The meaning of the subscripts on the variation should be clear, and will be dropped 
in the following. If we decide to constrain the model such that the trace of the matrices is to remain fixed (i.e. the 
variations preserve the trace of the matrices- a conclusion that as we shall see further on, will be forced on us 
through the spontaneous breaking of translational symmetry in the infinite N limit) then our variations must also 
have a vanishing trace. For the case where we consider variations of the form (|39|l (we shall consider linear 
combinations of such variations shortly) , this implies: 



= 



fdqSX,{q,q) = (e + e*)SaSiu - v) (40) 



This condition is automatically satisfied for variations of off-diagonal matrix elements. However for variations of 
diagonal elements, the phase factors would have to be purely imaginary for this expression to vanish. This will prove 
to have striking consequences further on. Note that these considerations would also apply to variations of finite 
dimensional matrices, however they do not alter the problem in quite the same way as we will see later. Indeed the 
"central extension" to the model alluded to earlier is a phenomenon unique to working with infinite dimensional 
matrices. 

At this point, it seems reasonable to question when, if ever the trace preserving constraint on our model is 
applicable since we claim such drastic consequences to follow from it. It turns out that this condition is redundant 
for finite dimensional matrices as our equations of motion view the trace as a center of mass co-ordinate for the 
collection of DO branes which moves with a constant velocity, and hence can be set to zero without loss of generality. 
This observation arises from taking the trace of the equations of motion- the trace of Xi gives us the acceleration of 
the center of mass, whereas the trace of the right hand side of the equation, involving only commutators vanishes. 
This is not the case in infnite dimensions. 

Now the trace of an unbounded operator is a slippery concept. In general it is going to depend on either a 
rcgularization scheme or an ordering prescription (equivalent to the concept of principal value in Riemann 
integration). To see this, consider the operator r whose spectrum as we have seen in the integers. The trace of this 
operator is only zero if we count "outwards from zero" . That is, 0-1-1 — 1-1-2 — 2-1-3 — 3... = 0. If we were to 
evaluate this same sum except now "counting outwards" from any other integer, we'd get an infinite number. Thus 
it would seem to be prudent at the very least, to restrict our variations to preserve the trace of our operators, 
whatever they may be, as once we've picked a prescription in which to make sense of the trace, arbitrary variations 



of the matrix elements are likely to derail this prescription. 

However, this argument is not very compelling from a physical perspective, and this need not be the case. 
Intuitively, one would expect that in the limit of a truly infinite number of DO branes, the center of mass co-ordinate 
will cease to be dynamical in the center of mass frame. The action |^ is invariant under a constant translation of 
the center of mass coordinate Xi —fXi + Cil, hence it posseses the 3-d translation group as a symmetry. If this 
symmetry is spontaneously broken in the limit of infinite dimensional matrices, then we can conclude that the center 
of mass co-ordinate (hence its velocity, as we chose to work in the centre of mass frame) is kinematically 
superselected to some fixed value and is hence a dynamical invariant. We thus have a new constraint on our problem 
that we had better respect when applying the variational principle- that our variations preserve the overall trace of 
the operators. We explore this possibility in the appendix, in what is essentially an application of the technique of 
the Coleman-Mermin- Wagner Theorem to our particular model [ig [l3 UM ■ 

However we will also see further on that the central extension we are about to derive, can also be derived from 
another fundamental requirement- that of Galilean invariance. This invariance arises in the limiting infinite 
momentum frame in which this matrix model is formulated [l9j. One finds that requiring this symmetry to be 
preserved in the limit of infinite dimensional matrices also induces a central extension to the model. 

Armed with this knowledge, we can proceed to vary the constrained action. The variation of Q by an arbitrary 
variation is given by: 

SS^ To J dtdqdp{^dX,{q,p)X,{p,q) + ^X,{q,p)SMp,q) ~je,jkSX,{q,p){[X,,Xk]){p,q)} (41) 

+ To J dtdqdpdr{^SX,iq,r)X,{r,p){[X,,X,])ip,q) - ^X,iq,r)SX,ir,p){[X,,X,]){p,q) (42) 

+ ^{[X,,X,]){q,p)SX,{p,r)Xj{r,q) - ^{[X,, X,]){q,p)Xj{p,r)dX,{r,q) (43) 

- -Y'^vkX^i<l,P)^XJiP^r)Xkir,q) + —e,jkX^{q,p)Xkip,r)6Xj{r,q)j (44) 

And the variation of the constraint term is given by: 

dt\{ f dqdpSX,{q,p){[[X,, [X„X,]]){p,q) - SX,{q,p){[X,, [X„ X,]]){p,q)} (45) 



SS = - 

By inserting the explicit form for the variations (|39|l . and taking care to respect the order of the integrations, 
extremizing the combined action with the Gauss law constraint term with respect to the variation 5i^^^u,v implies the 
following: 







i[-Xi + {[X,,[XuX,]])^znei,k{[X,,Xk])^2\{[Xu[X,,X,]])-2j^{[Xu[X,,X^^^^^^ (46) 



Where the second expression in the curly brackets differs from the first only through an interchange in the variables. 
As the equations of motion for the auxiliary field A imply that the last two terms in the above expression vanish, we 
can drop them right away. Recall that when u ^ v, the phase factor e is arbitrary. Hence the above must hold for 
e — \ and e = i, which is enough to guarantee the vanishing of each of the two expressions contained in the larger 
brackets separately, from which we conclude: 

Xi{u,v) ^ {[Xj,[Xi,X,]]){u,v) - ineijk{[Xj,Xk]){u,v) ; u^v (47) 

However when u = v, were our variations to be restricted by the trace preserving condition H40() . we see that the 
phase factor for this variation must be purely imaginary. We can see that in this case, (|46|l vanishes identically and 
so the principle of least action does not tell us anything. The most we can then conclude as our equations of motion 
is the following: 

Xi{u, v) = {[X,., [XuXj]]) [u, v) - iKeijk{[Xj,Xk])iu, v) + fi{u)5{u - v) (48) 



Where the fi{u) are some as of yet undetermined functions. Translating this expression back into operator language, 
we have the equations of motions: 

Xi - [X, , [Xi ,Xj]]~ ineijk [X, , Xk] + A; (49) 

Where the A; are diagonal operators. It is a straight forward exercise to show that we further constrain the A; to be 
multiples of the identity, upon considering more general linear combinations of H39|l that satisfy the trace preserving 
constraint. However a quick and hueristic way of seeing this is to consider an important aspect of the action Q that 
we've neglected up to now- it's U(N) global symmetry. In fact this virtue of the model alone, has generated much 
theoretical interest through the observation that as N goes to infnity, the symmetry group U{oo) could tend to the 
group of diffeomorphisms of certain 2-dimensional surfaces (see JJj and references therein) . This association 
appears in context of regularizing membrane theory by recasting it as a matrix model, where the original worldsheet 
diffeomorphism invariance of the theory translates into a U{N) symmetry |15J . More intriguingly, there are 
indications jlJl that C/(oo) is much bigger than this, and could contain as subgroups the different diffeomorphism 
groups of toplogically distinct manifolds, suggesting that a matrix model with this symmetry could be a theory of 
dynamical D-branc topology change. We investigatite this possibility hands on in Q (see |2j| for a similarly spirited 
investigation). 

Returning to the problem at hand, we can exploit this U{oo) invariance to deduce that the operators A; have to be 
multiples of the identity. 

A/ = cil (50) 

We stress that this particular argument is somewhat hueristic, but the result is nevertheless true from considering 
arbitrary traceless variations of the action (as dicussed above). Now in the case where N is finite, all one has to do 
is to take the trace of (|49|l to see that the c; all vanish. However, the fact that the trace of a commutator famously 
does not vanish when we are dealing with infinite dimensional matrices means that the c/ are now determined by the 
fact that overall, the trace of the right hand side of (|49|l has to vanish by virtue of the fact that the trace of Xi has 
to vanish. It is this fact that will permit us to look for new solutions to this matrix model in the limit of infnite 
dimensional matrices. For starters, we can take the c; to parametrize new solutions which may be central extensions 
of the usual algebras taken to define non-commutative geometries. However the most obvious new solution 
permitted is the non commuting plane- the equations of motion (|10|lHll|lH12|l now become: 

X+ = ^[X+,[X+,X^]] + [X3,[X+,Xs] + 2kX+] + c+I (51) 

X- = ^[X_,[X_,X+]] + [X3,[X_,X3]-2kX_] + c_/ (52) 

I3 = i[X+,[X3,X_]] + i[X_,[X3,X+]] + K[X+,X^]+C3l (53) 

From which we see that in the static case, for c± = 0, C3 = /?: 

X3 = A/ , X+^a'f , X_ = a ; [a, a^] = - (54) 

also solves our equations of motion (the case where all the Ci are zero would give us the fuzzy cylinder and the fuzzy 
sphere solutions). Here, A is a constant describing the displacement of the fuzzy plane along the z-direction. Hence 
in the limit of an infnite collection of DO branes in a RR background, a D2 branc forms with the structure of the 
non-commutative plane in addition to the fuzzy cylinder and fuzzy two sphere solutions. The non-commutativity of 
course being parametrized by the ratio of C3 over the magnitude of the RR field strength given by k. Interestingly, 
unlike the cases hitherto studied, the degree of non-commutativity dependends inversely on the background 4-form 
field strength. It should be clear that the Ci now serve to parametrize different matrix models, quite like the central 
charge of a given conformal field theory. We can of course pick certain values for them and look for new solutions. 



10 

but it would be desirable to understand the q in terms of more fundamental considerations. Such an understanding 
is readily available and we investigate this in the next section. However, to round out this part of the discussion, we 
study the energetics of these infinite dimensional solutions. It is straightforward matter to find the energy density e 
( = E/Tr I) of the fuzzy cyhnder: 

2 
^cyi = To-k'^p'^ (55) 



And the energy density of the fuzzy plane: 






^pi = ^0^75 (56) 



Where we should point out that the form of the Hamiltonian is no different in the infinite dimensional case. That 
the central extensions do not alter the energy of our solutions will be justified in the next section. As these are 
positive energies, one cannot interpret these solutions as bound states just yet, as one would have to study the local 
energy landscape around these solutions to test their (meta) stability. Such a study would have taken us far afield 
and we leave this to a follow up report. Since this concludes the main thrust of this report, we feel that a brief 
summary of our treatment is in order at this time. 

In studying the matrix model defined by (^ , we realised through naive manipulation of the static solutions for finit 
N that in addition to admitting fuzzy 2-sphere geometries as fundamental solutions, it appeared that fuzzy cylinders 
could also be solutions. Such solutions would neccesarily have been infinite dimensional as the lie algebra defining 
the fuzzy cylinder is that of the euclidean group in 2-dimensions, and as such is a non semi-simple algebra. We 
carefully considered the transition to infnite dimensional matrices in our model, and not only confirmed that the 
cylinder was indeed a solution, but discovered the fact that in general this model generically has central extensions. 
These central extensions are facilitated through the spontaneous symmetry breaking of translational invariance of 
this model in the limit of infinite dimensional matrices, and as we are about to see, through requiring that Galilean 
invariance be preserved in this limit. These extensions parametrize different models, and for special cases we are 
likely to find new fuzzy geometries as solutions. We considered a special central extension that permitted the 
non-commuting plane as a solution, thus further confirming the richness of physics contained in this matrix model. 
We hope to use the results here in a future report to investigate the rather pregnant possibility of dynamical 
topology change P| that this model seems to suggest as possible. 

IV. PHYSICAL INTERPRETATION OF THE CENTRAL EXTENSION 

We show in the appendix how the central extensions just uncovered, owe their origins to the spontaneous breaking 
of translational symmetry of the center of mass degree of freedom in the limit of infinite dimensional matrices. It 
should be clear that these central extensions can also be obtained by imposing the constraint that the center of mass 
degree of freedom be fixed through the method of Lagrange multipliers. Such a constraint can be enforced through 
the addition of the term: 

Or.nnatr = I Ot <7i\ — — - 1 rXi — A- 



j dt^Tra,{X,~\J) 



where the Ui are the functions of time which act as Lagrange multipliers, and the Xi describe values that the c.o.m 
coordinates are constrained to. Such an addition to the action will add the following "central" term to our equations 
of motion: 

X, = [X,,[X,,Xj]]-ie,jk{Xj,Xk]+cJ ■ c, = ajTrI (57) 

Where we see that the Ci we derived earlier are none other than the appropriate Lagrange multipliers rescaled. We 
expect the Ci to be non zero despite the factor of TrI in the definition, as the Lagrange multipliers are also expected 
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to scale proportional to TrI. This observation arises from the interpretation of the Lagrange multipUers as the 
appropriate force term enforcing the constraint, which we expect to scale with the number of degrees of freedom in 
the problem [20| ■ We should note here that this understanding of the Ci allows us to conclude that it would not 
contribute to the Hamiltonian of our model, as its constribution would vanish through the equations of motion for 
the lagrange multipliers. We should also note that inspitc of trying from the outset to avoid introducing unncccesary 
auxiliary fields, we see that for this particular constraint, it could not be avoided. This realisation arises when 
searching for the physical interpretation of the q. 

Having made this connection, one might also wonder what the effect would be of imposing as a constraint, 
something that one might ordinarily take forgranted from the Galilean invariance inherent in the model. Starting 
from the begining, let us add instead, for example the following term to our action: 

= I dt (Ji(——TrXi - u, 

dt ——-Trai{Xi - loJ) 

Such a term implies that the center of mass momentum is a constant. One would expect such an introduction to be 
redundant, by Galilean invariance. However examining the equations of motion without this term introduced into 
the action: 

X, = [Xj, [X,, X,]] - ie,,k[Xj,Xk] (58) 

We see that in the case of finite dimensional matrices, Galilean invariance is trivially implied for the center of mass 
momentum. However since it cannot be taken forgranted that commutators have a vanishing trace in infinite 
dimensions, we see that one might have a violation of Galilean invariance unless a compensating term appears in the 
equations of motion to account for this peculiar property of infinite dimensional matrices. Such a term is readily 
provided by the constraint term introduced above: 

X, = [X, , [X, , X,]] - ie.jk [X, , Xfe] +cj : c, = -^,/TrI (59) 

Where the only difference with H57|l is in the definition of the Ci . Thus we see that the central extensions to this 
matrix model can be deduced from seemingly disparate, but equally well grounded physical considerations. 

V. MORE SOLUTIONS 

In addition to the fuzzy plane solutions H54() . we find that the same central extension (parametrized by c± = 0, 
C3 = /3): 

X+ = ^[X+,[X+,X^]] + [Xs,[X+,X3]+2kX+] 
X- = i[X_,[X_,X+]] + [X3,[X-,X3]-2kX^] 
X3 = ^[X+,[X3,X^]] + i[X_,[X3,X+]] + K[X+,X^]+pi 



is also solved by the algebra: 



[X3,X±] = ±2kX± , [X+,X_] = ^I (60) 

K 



One can easily furnish a representation of this algebra through the identifications: 



X+ = a , X_ = a^ ; [a, a^] = - (61) 

K 

X3 = ^atfl (62) 
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Where we have reversed the identifications in (|54|l with X_|_ — a and X- = a^ . Although from H60() one might view 
this as a deformation of the algebra of the cylinder, it has a clearer geometric interpretation as a deformation of the 
fuzzy plane H54|l found earlier. One obvious reason for this is that the usual representation of the oscillator algebra is 
spanned by state vectors indexed by the non-negative integers, whereas the cylinder algebra is spanned by state 
vectors indexed by all the integers [26j- Furthermore (|60l) can be interpreted easily in terms of the fuzzy plane, 
where instead of being restricted to a surface of constant X3, it parabolically "warps" the further "out" you go 
(recall that X^X^ can be interpreted as a radial co-ordinate, where the identification is exact in the limit of 
commuting X^ and X_). Since X3 is quantized in integer steps of P/k, it would seem that this solution has warped 
the fuzzy plane into a parabola-like fuzzy configuration. 



VI. CLOSING REMARKS 

We conclude our report at this point, with the hope that the new solutions we have uncovered serve to demonstrate 
several (of the likely many possible) novel features of this model in the limit of infinite dimensional matrices. We 
feel that this may yet be a shadow of things to come, as previously in the discussion we raised the possibility that 
this model may be a candidate model of topological dynamics in the large N limit. We will investigate this 
possibility in a follow up report j^l j where transitions between the fuzzy cylinder and the (multiple) fuzzy sphere 
solutions will be investigated. Independent of this however, we feel that an equally fruitful avenue for future work 
would be to study the string theoretical interpretation of these results, which we have so far not attempted since we 
have deliberately chosen to study this model from the perspective of a non-commutative geometer. Certainly such a 
study is likely to uncover yet more of the richness of this model and perhaps shed further light on the some of 
questions that motivated this study, such as that of topological and geometrical dynamics in matrix models. 
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APPENDIX A: SPONTANEOUS SYMMETRY BREAKING AS N ^ 00 

Consider the action for the center of mass degrees of freedom, Xi — "^^' : 



S^To dt TrI x,x, (Al) 



Where we have used the relation that the kinetic energy = p^^^/2m. Now consider the quantum mechanical 
transition amplitude bewteen two states that differ by some center of mass displacement: 

(f(r)|f(0)) = f Dc e-^° ^0 '^^ ^^^^^"-^ (A2) 



Where we are functionally integrating over all histories that begin at a;(0) and end at ^(r) with d = x{T) — ^(0). 
Upon integration, any such amplitude will contain a prefactor coming from the classical path interpolating between 
these two states: 

x{t) = f (0) + x^ (A3) 
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Which will take the form (since TrI — N) 

{x{T)\x{0)) = e-^°^^{...) (A4) 

Where the remainder of the amplitude comes from integrating the fluctuations around this solution. From the 
prefactor, we can immediately see that this amplitude is vanishingly small in the limit iV — > cx). The only way this 
amplitude can have a finite prefactor is if the time taken for this interpolation, T becomes very large or if the 
displacement d becomes vanishingly small. This implies that the center of mass coordinate has to move 
infinitessimally slowly over very small distances. This is the same situation as if the center of mass were fixed. Hence 
we can conclude that the translation symmetry of ^ is spontaneously broken in the limit of infinite dimensional 
matrices. The displacement of the center of mass coordinate of our gas of DO branes is thus superselected to be 
some fixed value and in applying the variational principle, we should respect this superselection rule. 

It is instructive to make a comparison with the equivalent result in field theory (which illustrates the general 
conclusion that there are no Goldstone bosons in 1+1 dimensions [l^l)- Consider a free scalar field theory in 1+d 
dimensions: 

5'= i f dtd'^x d^,(j)d^(j) (A5) 

In the absence of a potential energy term, the value of the field serves as a modulus labelling different vacua. 
Consider an instanton transition between two different vacua given by the solution: 

</>(t,f)=(/.o + -^((/./-</)o) (A6) 

where tf is the time taken to complete the transition. Evaluating the action functional of the classical path 
interpolating between the two vacua gives us: 

S,^\^^^l^V (A7) 

Thus if we begin by studying this theory in a d-dimensional box of length L, and take the limit t^ — > oo, L — > oo 
isotropicaly, we see that the classical prefactor multiplying the quantum transition amplitude becomes: 



-\(4>f-<po)^L 



(A8) 



From which we see that if d > 2, this transition amplitude vanishes and we have spontaneous symmetry breaking of 
the translation symmetry in (j). However when d = 1, this amplitude remains finite and hence transitions between 
the different vacua are generically possible, hence this symmetry is not spontaneously broken (no goldstone bosons 
in 2-d). 
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